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Abstract . One of the main aims in the next generation of precision cosmology expe- 
riments will be an accurate determination of the equation of state (EOS) for the dark 
energy (DE). If the latter is dynamical, the resulting barotropic index uj should exhibit 
a non-trivial evolution with the redshift. Usually this is interpreted as a sign that the 
mechanism responsible for the DE is related to some dynamical scalar field, and in some 
cases this field may behave non-canonically (phantom field) . Present observations seem 
to favor an evolving DE with a potential phantom phase near our time. In the literature 
there is a plethora of dynamical models trying to describe this behavior. Here we show 
that the simplest option, namely a model with a variable cosmological term, A = A{t), 
leads in general to a non-trivial effective EOS, with index Wg, which may naturally 
account for these data features. We prove that in this case there is always a "crossing" 
of the We = — 1 barrier near our time. We also show how this effect is modulated (or 
even completely controlled) by a variable Newton's gravitational coupling G = G{t). 
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1. Introduction 

Cosmology is becoming a science of precision and accuracy. In the last few years a flurry of 
experimental activity has been devoted to the measurement of the cosmological parameters. This 
effort has transformed observational cosmology into a fairly respectable branch of experimental 
physics ^ m Ol- As a result, evidence is piling up in favor of the existence of the dark energy 
(DE) component, pD, pervading the Universe (or at least the known patch of it). Hopefully, this 
situation is bound to improve further with the advent (in the near future) of the promising SNAP 
and PLANCK projects which will grant an experimental determination of the cosmological 
parameters to an unprecedented few percent level of accuracy. Quite in contrast, the nature of 
the DE remains at the moment a profound mystery. Historically, the DE was first identified with 
the cosmological constant (CC), A, and the vacuum energy contributions it receives from quantum 
field theory (QFT) [S]. This idea led to the famous cosmological constant problem [SI EI. In more 
recent times the notion of DE has been detached from that of A and has been extended to a 
variety of models leading to an accelerated expansion of the universe in which the DE itself is a 
time-evolving entity. These models include dynamical scalar fields (quintessence and the like) [SI El j 
phantom fields |lUj , braneworld models , Chaplygin gas , and many other recent ideas like 
holographic dark energy, cosmic strings, domain walls etc (see e.g. [7] and references therein). 
Obviously the very notion of CC in such broader context becomes diminished, the CC could just 
be inexistent or simply relegated to the status of one among many other possible candidates. 

A general DE model is described as being a sort of fluid characterized by a conserved DE energy 
density, po, and pressure related by an equation of state (EOS) po = ojd Pd- A particular and 
popular realization of the DE is the aforementioned quintessence idea [SI El) where one has some 
scalar field x which generates a non- vanishing /9£) from the sum of its potential and kinetic energy 
term at the present time: p^ = {(1/2)^ + V{x)}t=to- The sign of the coefficient ^ determines 
whether the field can describe quintessence (^ > 0) or phantom DE (^ < 0). If the kinetic energy 
for X is small enough, it is clear that p^ looks as an effective cosmological constant. The scalar 
field X is in principle unrelated to the Higgs boson or any other field of the Standard Model (SM) 
of particle physics, including all of its known extensions (e.g. the supersymmetric generalizations 
of the SM); in other words, the x field is an entirely ad hoc construct just introduced to mimic the 
cosmological term. However, the field x is usually thought of as a high energy field (unrelated to 
SM physics), i.e. x — where Mx is some high energy scale typically around the Planck mass 
Mp ~ lO"*^^ GeV. If one assumes the simplest form for its potential, namely V{x) = (1/2) x^j 
its mass turns out to be m^^ ~ ifo ~ 10~^^ eV in order to approximately describe the present 
value of the energy density associated to A, namely = Ao/Svr G ~ 10"^^ GeV^. Obviously it is 
very difficult to understand the small mass m^^ in particle physics, and this is of course a serious 
problem underlying the quintessence models. On the other hand the possible observed transition 
from a quintessence into phantom-like behavior jl3l I14j cannot be explained with only one scalar 
field (because the sign of the coefficient ^ is fixed), and therefore one has to resort to at least two 
scalar fields, one with > and the other with another coefficient ^' < ^. In spite of the many 
^For various recent theoretical approaches to the lod = —1 boundary crossing see e.g. |15f . 
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difficulties with the scalar field models they have the virtue that they may help to understand why 
the DE might be evolving with time as suggested by the analyses of the cosmological data jl3l I14j , 
and in this sense one would like not to lose this useful aspect of the quintessence proposal. 

In this work, in contrast to the previous approaches, we wish to stick to the original idea that the 
primary cause of the DE is a "true" CC term in Einstein's equations, but we assume it is a variable 
one: A = A(t). Actually, a rich variety of variable A models have been reported in the literature - 
some of them admittedly on purely phenomenological grounds jl6l I17j . There are however models 
based on a more fundamental premise, even if ultimately phenomenological at the moment. For 
instance, there are models attempting to adjust the small value of the cosmological term with the 
help of a dynamical scalar field ? as for example the cosmon model [SI E| • More recently a 
variable A has been proposed as a running quantity from the point of view of the renormalization 
group (RG) |2Uj-|24j. Particularly interesting for the present work is the possibility that an effective 
equation of state can be associated to a variable cosmological model. This has been explored in 
detail in Ref . for a model of running cosmological term |^ . In the present work we generalize 
this idea to arbitrary variable cosmological models (without specifying the underlying dynamics), 
and show that these models can effectively lead to a non-trivial EOS (therefore mimicking a 
dynamical DE model). We also show how to include in the effective EOS the possible effects 
from a variable Newton's gravitational coupling. We provide a general algorithm to construct this 
effective EOS and prove some general results which describe the relationship between the variable 
A models and the dynamical DE models. Finally we argue that this relationship may help to shed 
some light to understand the apparently observed transition from quintessence-like to phantom-like 
behavior as indicated by some analyses of the most recent cosmological data |131 114j . 

2. Two cosmological pictures: Variable A and G versus dynamical dark energy 

We start from Einstein's equations in the presence of the cosmological constant term, 

^i^y - ^9tJ.uR = 8ttG T^v + g^lvI^ = 8ttG {T^y + g^^ Ph) , (1) 

where T^,y is the ordinary energy-momentum tensor associated to isotropic matter and radiation, 
and Pa represents the energy density associated to the CC. Let us next contemplate the possibility 
that G = G{t) and A = A(i) can be both functions of the cosmic time within the context of 
the FRW (Friedmann-Robertson- Walker) cosmology. The possibility of a variable gravitational 
coupling and a variable cosmological term in cosmology (variable CC term, for short) has been 
considered by many authors, and in particular it has been discussed within the renormalization 
group approach to cosmology |2flj-|24j. It should be clear that the very precise measurements of G 
existing in the literature refer only to distances within the solar system and astrophysical systems. 
In cosmology these scales are immersed into much larger scales (galaxies and clusters of galaxies) 
which are treated as point-like (and referred to as "fundamental observers", co-moving with the 
cosmic fluid j26j). Therefore, the variations of G at the cosmological level could only be seen at 
much larger distances where we have never had the possibility to make direct experiments. To put 
it in another way, the potential variation of G = G{t) and A = A(t) is usually tested in terms of a 
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measurable redshift dependence of these functions, G = G{z) and A = A(2;), say a variation for a 
redshift interval of at least Az > 0.01, therefore implying scales of several hundred Mpc. 

The only possible EOS for the CC term, whether strictly constant or variable, is p\ = —p\. 
However, we may still parametrize this variable A model as if it would be a dynamical field 
DE model with pressure and density {pd,Pd)- We will call this the "effective DE picture" of the 
fundamental A model. The matching of the DE picture with the fundamental CC picture generates 
a non-trivial "effective EOS" [Uj for the latter, pD = LUePn, where the effective barotropic index 
(Je = uje{z) is a function of the cosmological redshift to be determined. By exploring this alternative 
we can test the variable A models as a canonical source of effective dynamical DE models. In the 
following we will restrict our considerations to the flat space case {Qk = 0). This is not only for 
the sake of simplicity, but also because it is the most favored scenario at present. Friedmann's 
equation with non-vanishing p\ reads 

^1 = ^(p + Pa), (2) 

where for convenience we have appended a subscript in the Hubble parameter. On the other hand 
the general Bianchi identity of the Einstein tensor in leads to 

VMG'(T^. + 5m.Pa)] =0. (3) 

Using the FRW metric explicitly, the last equation results into the following "mixed" local conser- 
vation law: 

^ [G{p + Pj,)]+3GHa{p + p) = 0. (4) 

If PA = dpx/dt 7^ 0, /3 is not generally conserved as there may be transit of energy from matter- 
radiation into the variable p\ or vice versa (including a possible contribution from a variable G, if 
G ^ 0). Thus this law indeed mixes the matter-radiation energy density with the vacuum energy 
(pa)- To be more precise, the following scenarios are possible: i) G =const. and p\ =const. 
This is the standard case of ACDM cosmology implying the canonical conservation law of matter- 
radiation: p + 2>H [p + p) = 0, ii) G =const and p\ ^ 0, in which case Eq.(IlJ) boils down to 
PA + P + 'iH (p + p) = 0; iii) G ^ and p^ =const, implying G{p + pa) + G[p + 3H{p + p)] = 0; 
and finally iv) G 7^ and pA 7^ 0, which in the case of self-conservation of matter-radiation it 
leads to {p + pa)G + GpA = 0. Notice that in cases ii) and iii) the matter-radiation cannot be 
canonically conserved (if one of the two parameters pA or G indeed is to be variable), whereas in 
case iv) it is assumed to be conserved. Explicit cosmological models with variable parameters as 
in cases ii) and iv) have been constructed within the context of RG cosmology in Ref . |22j and j23j 
respectively. 

With only Eq. (^J and the FRW equations we cannot solve the cosmological model with variable 
parameters. One needs a fundamental model that informs us on the functional dependence of pA 
and G on the remaining cosmological functions, like p,p and H. However, in most of this work 
we will not commit ourselves to any specific model for the underlying fundamental dynamics, 
because the kind of results we are aiming at are to be valid for a large class of models with 
dynamical cosmological parameters. Our guiding paradigm will be the general expectation that 
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these parameters are variable because of the effective running that they should display according 
to the renormalization group. In this sense the possible variation of "fundamental constants" such 
as A and G could be an effective description of some deeper dynamics associated to QFT in curved 
space-time, or quantum gravity or even string theory, all of which share the powerful RG approach 
to quantum effects. For instance, in the specific framework of QFT in curved space-time the RG 
equations for p\ and G should entail definite laws pA = PK{t) , G = G{t) pSl. Given a fundamental 
model based on QFT, the parameter p\ will primarily depend on some cosmological functions 
(matter density p, Hubble expansion rate H, etc) which evolve with time or redshift. Similarly for 
the Newton coupling G. Therefore, in general we will have two functions of the redshift 

Pa{z) = pMz),Hiz), ...) , Giz) = Gipiz), Hiz), ...) . (5) 

We understand that other fundamental parameters could also be variable. For example, the 
fine structure constant has long been speculated as being potentially variable with the cosmic 
time (and therefore with the redshift) including some recent experimental evidences - see e.g. 
j29j . This might also have an interpretation in terms of the RG at the level of the cosmological 
evolution. And similarly with other "constants" such as the ratio nie/mp between the electron 
to the proton mass etc which could also have evolved throughout the cosmic history. However, 
in this paper we concentrate purely on the potential variability of the most genuine fundamental 
gravitational parameters of Einstein's equations such as A and G. The possibility that G could 
be associated to a variable scalar field stems from the old Jordan and Brans-Dicke proposals [HO] 
and has generated an abundant literature since then. Similarly the idea that the A term could 
perhaps be variable and related to a dynamical field is also relatively old |181 1191 In short, 
whether related to dynamical fields or to the general phenomenon of RG running, the so-called 
"fundamental constants of Nature" are suspicious of being non-constant from the point of view of 
quantum theory. Adopting this general Ansatz, we wish to investigate whether general (model- 
independent) properties can be extracted if we take Einstein's equations in the FRW metric, 
with variable parameters © , as the starting point for the study of cosmology. 

Functions © will usually be monotonous; e.g. in the RG model of Ref. ^27 p\ inherits its 
time/redshift dependence through pA = ci + C2 and it satisfies dp^/dz ^ if C2 ^ respec- 
tively. Similarly, G = Gq/{1 + i^lnH^/H^) in the RG model of and therefore dG/dz ^ 
if ^ respectively. Furthermore, in most phenomenological models in the literature (where 
the underlying dynamics is usually not specified) this monotonic character also applies P^-jl9j. 
Since the functions © are presumably known from the fundamental model, and the ordinary EOS 
for matter-radiation is also given, one may solve Eq. @ and Friedmann's equation to determine 
p = p{z), p\ = pa{z) and G = G{z) as explicit functions of the redshift. These may then be 
substituted back in to get the expansion rate of the variable A model as a known function of 
z: 

Hl{z)=Hl [Ql, fM{z;r){l + zT + nlfA{z;r)\ , (6) 

where a = 3(1 -|- Um) {^m = 1/3 or tOm = for the radiation or matter dominated epochs 
respectively). Functions Jm and /a are completely determined at this stage, and may depend on 
some free parameters r = ri, r2, .. - see [221 123| for non-trivial examples with a single parameter. 
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Only for the standard AC DM model we have /m = /a = 1- In general it is not so due to the 
mixed conservation law ((!]). Whatever it be their form, these functions must satisfy /M(0;r) = 
/a(0; r) = 1 in order that the cosmic sum rule + fi^ = 1 is fulfilled. 

On the other hand, within the context of the effective DE picture the matter-radiation density 
satisfies by definition the standard ("unmixed") conservation law 

Ps + aHDPs = 0. (7) 

Here we have denoted the matter-radiation density by ps to emphasize that in the DE picture 
it satisfies the standard self-conservation law, i.e. it corresponds to case i) mentioned above. 
Moreover Hd stands for 

Hh = '-^iPs + PD) (8) 
with constant Gq. The subscript D here is to distinguish it from the Hubble function in the CC 
picture, Eq.©. The DE density pn is also conserved, independently of matter: 

PD + 3HD{l+UJe)pD=0- (9) 
The solutions of the two conservation equations (O and ^ read ps{z) = Ps{0) (1 + z)" and 

PD(.z)=pD{0)az), (10) 

From it is clear that if 0Je{z) > —1, then po increases with z (hence decreases with the 
expansion). This is characteristic of standard quintessence |31j. In contrast, if at some point 
z = Zp we hit 0Je{zp) < —1, the DE density will decrease with z and hence will increase with the 
expansion. This anomaly would signal the breakthrough of a phantom regime ^H] in some interval 
of z around Zp. The general solution for ps{z) and pd{z) given above can be substituted back to 
Friedmann's equation in the DE picture and one gets the corresponding expansion rate: 

hUz) = Hi[nl,{i + zr + n%az)]. (11) 

Here + Cl^ = 1 because C(0) = 1. In general AQm = - ^li 

is nonzero. In fact, suppose 

that one fits the high-z supernovae data using a variable p\ model - cf. the second Ref . |22j . 
The fit crucially depends on the luminosity distance function, which is determined by the explicit 
structure of ©, so that the fitting parameters ri^j,r2^ can be different from those obtained by 
substituting the alternate function (jllj) in the luminosity distance. 

3. Effective EOS for a general model with variable A and G 

Let us now compute Ue- With the help of ([T?H) we have uje{z) = — 1 + (1/3) ((1 + z)/() dQ/dz. 
Differentiating Eq. Hll|) on both sides we can eliminate dQ/dz. Next we apply the matching con- 
dition mentioned above, meaning that we fully identify with H\ in the effective DE picture. 
The final result is 

-l + {{l + z)/^){l/Hl{z))dHl{z)/dz + {l-a/^)^l, {Hl/Hliz)) (1 + z)" 

^'^'^ i-^%{Hi/Hi[z)){i + zr ■ ^ ^ 
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This is the effective barotropic index for a variable pj\,{z) and/or G{z) model. The procedure is 
formally similar to Ref |32j but it is conceptually different. Here we do not try to reconstruct the 
scalar field dynamics as a fundamental model of the DE, not even a model-independent polynomial 
approximation to the DE jj^; rather, we start from a general (purportedly fundamental) rnodel 
that can decay into matter and radiation and then we simulate it with a (generic) dynamical DE 
model in which the matter-radiation is strictly conserved, as usually assumed. The rationale 
for this is that fundamental physics (e.g. the RG in QFT) can offer us useful information on 
the possible functional forms ©, as shown in[^-|21]. As a consequence, H\{z) can actually 
be computed after solving the coupled system of equations formed by Friedmann's equation, the 
non-trivial conservation law Q and the relations © provided by the fundamental model. Then 
one can check whether the effective EOS for the variable /oa model does emulate quintessence or 
phantom energy depending on whether the calculation of the r. h. s. of H12|) yields — 1 < < — 1/3 
or We < —1 respectively. 

It is remarkable that the tracking of the p\{t) model by the DE picture (in particular a possible 
phantom behavior near our time) is actually expected for a wider class of variable pA models. It 
is not easy to see this from Ea. (|12() . Instead, let us return to Eq.@. From the aforementioned 
matching condition we have G{p + pa) = Go{ps + Pd)- Using also H dt = —dz/{l + z) we can 
transform the general Bianchi identity (jlj) into the following differential form 



(1 + z) d{ps + pd) = a {ps + PD- ^a) dz , 



(13) 

where we have defined 

Uz) = ^pUz). (14) 

This new variable coincides with p\ only if G has the constant value Gq, but in general both p\ 
and G will develop a time/redshift evolution from the fundamental theory (for instance from given 
RG equations), so that will be a known function ^a = Chip, H, ...). Using ((TJ, rewritten in 
terms of the variable z, we can eliminate ps from H13|) and we are left with a simple differential 
equation for p£): 

dppjz) Pd{z) -U{z) _ 

— = a — = P{pd{z)). (15) 

dz I + z 

This equation is particularly convenient to analyze the behavior of the effective DE picture. Con- 
sider the plane {po, P{pd))- For constant ^a (as in the ACDM model, where G = Gq and = /^a)' 
Eq. H15|) resembles a renormalization group equation for the "coupling" po, with function" 
P{pd) and a fixed point (FP) at pj-, = ^Ai where (3{pJ^) = 0. Since dj3{pD) / dpo = a/{l + z) > 
(for all z > —1) it is clear that this FP is an infrared (IR) stable fixed point, hence Pd{z) p*d 
for decreasing z. For variable .^a = Ca(-z) the conclusion is the same: the solution pD = pd{z) of 
(|15j) is attracted to the given function ^a = Ca(-2) in the far IR. The quintessence regime occurs 
for (3{pd) > 0, i.e. for po > ^a, whereas the phantom regime is characterized by (3{pd) < 0. 
Whatever it be the sign of (3{pd) at a particular instant of time in the history of the Universe, 
the function pd{z) will eventually be driven to iK{z) for decreasing z (namely by the expansion of 
the Universe). Let us next study if the two functions Pd{z) and ^^{z) can be very close at some 
point near our time, and what are the phenomenological implications. The differential equation 
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(|15|) can be integrated in closed form for any (,Aiz)- 

Pd{z) = {1 + z) 

Expanding (,a{z) around z = —1 one can check that Pd{z) — > ihiz) at sufficiently late time (i.e. 
z — > —1). Finally, from (jlUp and with the help of (|15() we find 

^e{z) = -1 + f f 1 - ^) - -1 + , (17) 

3 V Pd{z)J 

where pd given by (|16() . This formula is more useful than Ea.()12|) in the present context. 
It provides an efficient recipe to compute We directly from the sole knowledge of ^a- Using these 
general formulae we can e.g. reproduce the effective EOS obtained for the particular model recently 
studied in Ref. j25j . see Eq. (23-24) of the latter. At the present time (that is, near z = 0) we expect 
that e{z) is a relatively small quantity because observations show that the effective barotropic index 
must be near —1 at present 0|2llS|. The sign of e{z) at z = 0, however, is not known with certainty. 
Model independent analyses of the most recent cosmological data seem to indicate |131 114j that 
cUe lias undergone a certain evolution from tJe < (at z ~ 1.7) reaching the phantom regime 
We ^ — 1 at z > 0. In j25j it was shown that a RG model for the p\ evolution can exhibit this 
type of behavior. We will now prove that this feature actually holds for a large class of models 
with variable .^a- 



Pd{0) 



a 



(16) 



4. Effective quintessence/phantom behavior of models with variable A and G 

Let us rewrite the solution of the differential equation (|15|) in the following alternative way: 

where z* is a root of I3{pd{z)) = 0, i.e. a point where pd{z*) = ^a(-z*)- Remarkably, one can show 
that a root z* always exists near our present time, meaning in our recent past, in the immediate 
future or just at z* = 0. The proof is based on establishing the following relation: 

dC{z) a(l + z)°-i 



dz 1 - nl, 



{nl,fMiz;r)-Ql) . (19) 



It ensues after a straightforward calculation from: i) the matching condition of the two pictures, 
H\{z) = H£){z), and ii) the constraint imposed on the functions Jm and /a in © by the Bianchi 
identity @. Since /M(0;r) = 1, as noted previously, it is clear that if the cosmological parameters 
of the two pictures coincide (AQm = 0) the derivative (|lljj) vanishes identically at z* = 0, hence 
a'e(O) = —1. If, however, AQm but is small (after all Cl^j and il^^ in the two pictures should 
not be very different), then there exists a point ^ not very far from z = where d(^/dz\z=z* = 
(and so uJe{z*) = —1). This completes the proof. Therefore, z* ~ exists and defines a (local) 
divide between a quintessence phase and a phantom phase. The slope of the function pD reads 

dppjz) a-i r dz'dijj^ 

; = — q(1 + z) / — — . (20) 

dz ^ ' (1 + z'Y dz' ^ ' 
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One could naively think that for increasing/decreasing with redshift, Ue{z) should always be 
above/below —1. Let us assume that is monotonous (hence dS,\/dz has a definite sign) in some 
interval containing z* and some reference point z = zi. Eq. ()2U() shows that in this case dpo/dz 
is also monotonous, and if z* < z\ then dpj:)ldz\z=z^ has opposite sign to d^^t^j dz\z=ZY- Thus e.g. 
if d^h/dz > (i.e. for decreasing with the expansion) then dpD/dz\z=zi < 0, meaning that 
the dynamical DE picture will look as phantom energy at zi (contrary to naive expectations); if 
z* > z\ then dpD/dz\z=zi > and the DE will look as quintessence at zi. Similarly, if S,a{z) is 
monotonically increasing with the expansion [d^p^/dz < 0), and z* < zi, the observer at zi will 
see quintessence (again counterintuitive), but when zi < z* he/she will see phantom DE. The last 
situation is particularly worth noticing because if = this case could just correspond to the 
present observational data (if the tilt into the phantom regime is finally confirmed |131 114j ) . The 
exact position of z* depends on the peculiarities of the variable ^a = (,a{z) model and also on the 
values of the parameter difference AQm in the two pictures © and This is corroborated 

in concrete scenarios with monotonically decreasing pa{z) _25j, where for parameter differences of 
a few percent one finds that the effective EOS develops a phantom phase at z* > and persists 
phantom- like asymptotically until the remote future, where ijOe{z) — > —1. 

Finally, we note from (jlj that if G = G{t) and we assume that matter is conserved (no transfer 
of energy with a variable pa), then dS,\/dt = —{p/Go) dG/dt. In this case if G{t) is monotonous, 
^A(i) is too. From (^0)) we get 

^ = a(l + z)"-i^[G(z)-G(z*)]. (21) 

This shows that if G is asymptotically free - that is to say, if G decreases with redshift - we should 
observe quintessence behavior for ^ z ^ z*; whereas if G is "IR free" (i.e. it increases with z - 
hence decreases with the expansion), then we should have phantom behavior in that interval. In 
general, if pA is variable (and G is fixed or variable) the monotonous property of ^a is satisfied 
by most models in the literature , including the RG ones |2n[ I2H I22[ I23| I24j . The monotonous 
variation of (,a{z) around z* is important to insure a long quintessence-like regime preceding this 
point - as suggested by the cosmological data. 



5. An example: renormalization group model with running p\ 

As a concrete example of cosmological model with variable parameters, consider the case G = 
const, and pA 7^ 0. We have catalogued this possibility in Sect. 2 as case ii), for which the general 
Bianchi identity (0)) simplifies into 

PA + p + 3H{p + p) =0. (22) 

This equation, in combination to Friedmann's equation, has been used to solve the cosmological 
RG model of j^. This model treats the CC semiclassically as a running parameter, therefore one 
evolving with time/redshift due to the quantum loop effects of the high energy fields (the only 
ones that can contribute significantly in this model, due to the "soft-decoupling" phenomenon 
associated to the CC in Ref. [^). Interestingly enough this model can simulate an apparent 
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phantom behavior near our time|2oj, actually following very closely the polynomial data fits of 
Ref. ^H]- In this particular model equations © read 

p^{z) = Co + CiH'^{z), G = const, (23) 

with 

Co = PA,o - ^Mj. Hi , Ci = ^Ml. (24) 

This model has a single parameter, z^, defined essentially as the ratio (squared) of the masses of 
the high energy fields to the Planck mass (Mp) j22j . 

a M2 



127r Ml 



(25) 



Here o" = ±1 depending on whether bosons or fermions dominate in their loop contributions to 
the running of p\. The typical value for u is obtained when M = Mp, namely 

vo = -^c^ 0.026 . (26) 

In general we expect \v\ ^ because from the effective field theory point of view we assume 
M ^ Mp. This is also suggested from the bounds on u obtained from nucleosynthesis and also 
from the CMB [SSI- 

Using equations (|22() and ()23() in combination with Priedmann's equation one can solve this 
cosmological model explicitly and one can obtain the functions H = H{z;v), p = p{z;v) and 
PK = Pk{z] v) in close analytic form. We refer the reader to [22] for all the details concerning this 
running CC model. Once the model is solved we have an explicit expression for the function ^a(-2^); 
Ea. ((Tl)) . which in this particular case is obviously given by ih{z) = pa{z; v). For the flat case one 
finds (in the matter-dominated epoch, a = 3) 

U{z; v) = pI + pI, \{\ + zt^-'^ - ij , (27) 

\ — V V J 

where and are the values of these parameters at z = 0. Notice that for = we recover 
the standard FRW case with constant cosmological term. Substituting ()27() in the general formula 
H16|) . and integrating, one obtains the effective DE density associated to this model: 

p^(z; v)=pU (pI - pI) (1 + zf + -^{l + zf [(1 + z)-3- - 1] + [(1 + zf - l] . (28) 

1 — V 1 — V 

As expected from the general discussion below Ea. HlSf) we confirm that poiz] v) — > iK{z] v) in the 
remote future; in other words uJe{z — > —1) — > —1, see Eq. H17|) . From (^7)) and (P5|l we have the 
asymptotic limit 

PD{z]iy)^U{z;v)^ pI-YZ^ (forz^-1). (29) 

Finally, the effective EOS parameter for this model is obtained by inserting (|27|) and (|28j) in 
Ea. ()17() . Upon some rearrangement it finally yields 

^ 170, [(1 + Z)3(l-) - 1] - (1 - ^) [J^O, (1 + Z)3 - 1] ' 
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(a) (b) 

Figure 1: (a) Numerical analysis of the effective EOS parameter Wg, Eq. (|17|) . as a function of the 
redshift for fixed v = —O.Suo < and for AQm = —0.01. The Universe is assumed to be spatially 
flat [Q^ = 0) with the standard parameter choice = 0.3 , = 0.7. In the figure, z* is the 
crossing point of the "barrier" uJe = —1; (b) The corresponding evolution and crossing (at z = z*) 
of the density functions S,a{z) and Pd{z) in the two pictures, Eqs. (|7rj) - H28() . in units of the critical 
density pcfi at present. 

where we have defined = p%/ Pcfi and Cl\^ = p^/pcfl corresponding to the matter densities in 
the two pictures. Recall that in the CC picture pM is non- conserved in this model - see Ea. iP^ 
- whereas in the DE picture ps is conserved - cf Eq Q -, so the two mass density parameters p^^ 
and /9g involved in the two fits (in each picture) of the same data need not to coincide. In our 
general discussion below Eq. (|11() we expressed this fact by asserting that the parameter difference 
Ail A/ = Q^.j — Cl^.j is in general expected to be nonzero^. 

In Fig.^ we show the numerical analysis of Eq. (|30jl for a typical choice of the parameters. 
In general we do not expect AQ]\f to be very large because the two pictures are supposed to give 
a similar representation of the same data. Therefore we have chosen AQm of order one percent. 
Similarly, the parameter should be close to 1)26^ . For A^Im < and < as given in Fig^ we 
find a transition point from quintessence-like behavior into phantom-like behavior very near our 
recent past, namely around z* ~ 0.7. - marked explicitly in the figure - i.e. around the time when 
the Universe deceleration changed into acceleration. In Fig.^ we plot the density functions H27|l 
and H28|) normalized to the critical density at present, pc,o- This plot also exhibits the crossing point 
z = z* oi the two curves, where Ca(-z*) = Pd{z*)- As predicted by the general differential equation 
((T3|) . in this figure the quintessence-like behavior is seen to be characterized by Pd{z) > S,a{z) 
whereas the phantom-like behavior by Pd{z) < Ca(-z)) with z = z* acting as a transition point 
between the two. We can also check in Fig.^ the accomplishment of the asymptotic condition 
given by Eq. ((^ . 

The interesting features displayed by this running cosmological model are only a particular 

illustration of the general expectations for cosmological models with variable cosmological param- 

■^We note that the result lljUI I. obtained from the general method devised in this paper, does perfectly agree with 
the particular calculation performed in I25| . 
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eters. The alternate representation of the effective DE density Pd{z) and its derivative dp^/dz - 
see equations ((T^ and (PU)) - allowed us to describe these features in a completely general way 
without committing to any particular model. If we change the model (namely the kind of evolution 
of the redshift functions ©) similar features will arise depending on the parameters of the model 
and on the value of Afijv/, but in all cases a transition point z* near our present time (lying either 
in our recent past or in our immediate future) will be found. The existence of the crossing point 
z* {uj(,{z*) = — 1) where the effective EOS may change from quintessence-like to phantom-like be- 
havior near our present time is, as we have proven in this paper, a general result for all the models 
with variable cosmological parameters. 

6. Conclusions 

To summarize, we have shown that a model with variable p\ and/or G generally leads to 
an effective EOS with a non-trivial barotropic index, with the property iVe{z) — s- — 1 in the far 
IR. Such model can mimic a dynamical DE model and effectively appear as quintessence and 
even as phantom energy. The eventual determination of an empirical EOS for the DE in the 
next generation of precision cosmology experiments should keep in mind this possibility. The 
nature of the dynamics behind /3a is not known in principle, but various works in the literature 
suggest that a fundamental px can display a RG running which can be translated into redshift 
evolution. Here we have generalized these ideas and have shown that in any model with variable 
^a(-z) = {G{z)/Gq) p\{z) there is an effective DE density, poiz), that tracks the evolution of (,a{z) 
and converges to it in the far IR. Moreover, we have proven that there always exists a point z* 
near z = where Pd{z*) = ^/^(z*), hence tOeiz*) = —1. If this point lies in our recent past {z* > 0) 
and .^A is a decreasing function of z around z*, then there must necessarily be a recent transition 
into an (effective) phantom regime LOeiz) < —1. This would explain in a natural way another 
"cosmic coincidence": why the crossing of = — 1 just occurs near our present time? Our results 
are model-independent, they only assume the formal structure of Einstein's equations of General 
Relativity and FRW cosmology in the presence of a variable cosmological term, p\ = p\{z), and 
possibly (though not necessarily) a variable gravitational coupling G = G{z). We conclude that 
there is a large class of variable {p\, G) models that could account for the observed evolution of the 
DE, without need of invoking any combination of fundamental quintessence and phantom fields. 
Last, but not least, we wish to point out that PLANCK and SNAP 0] data can constrain the 
underlying model with varying parameters (e.g. the RG model of Sect. 5) - see also j22| . Since, in 
principle, the underlying model may also have implications at astrophysical scales J^,, one could 
even check if the constraints from SNAP and PLANCK can yield the expected behavior at these 
scales. Most important, if the SNAP and PLANCK data should confirm the crossing of the CC 
boundary, u)e{z*) = —1, it would strengthen this approach substantially since the CC boundary 
crossing is generic in the entire class of cosmological models with variable cosmological parameters 
studied in this paper. Notice, however, that if observations would prove that Ue > —1 this would 
not invalidate the model because the crossing point could be in our immediate future {z* < 0). 
What we have proven, indeed, is that in this kind of models there always exists a crossing point 
z* around our present (z = 0), whether in our recent past or in our near future. 
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